The Spin Glass Transition : Exponents and Dynamics 
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Numerical simulations on Ising Spin Glasses show that spin glass transitions do not obey the usual 
universality rules which hold at canonical second order transitions. On the other hand the dynamics 
at the approach to the transition appear to take up a universal form for all spin glasses. The 
implications for the fundamental physics of transitions in complex systems is addressed. 



Spin glasses are magnetic systems with conflicting and 
random interactions between individual spins. The de- 
nomination "spin glass" for the wide class of such mate- 
rials which show a complex form of magnetic order at low 
temperatures dates from 1970 the same year and in- 
dependently it was found that a dilute magnetic alloy has 
a sharp cusp in the low field a.c. susceptibility at a tem- 
perature T g Q , so defining a precise critical temperature 
despite the disorder. Spin glasses have been intensively 
studied ever since, theoretically, experimentally and nu- 
merically in the belief that they will provide a paradigm 
for a vast family of complex systems. 

Numerical work has concentrated on the Edwards- 
Anderson version of a spin glass - a regular lattice of 
Ising spins coupled by random interactions. It turns out 
that this "theoretician's spin glass" shows the same es- 
sential characteristics as the experimental systems, and 
is a perfect test-bed for checking out ideas on complex 
systems in general. The mean field (or infinite dimen- 
sion) limit of this model is well understood However, 
in finite dimensions even this simplified version of a com- 
plex system turns out to be very devious, so that after 25 
years of work there is still little concensus on the funda- 
mental empirical facts, let alone the underlying theory. 

At low temperatures there is a domain of irreversible, 
non-ergodic behaviour; we will discuss exclusively the dy- 
namics in the other, ergodic, regime from high temper- 
atures down to (or nearly down to) the freezing tem- 
perature T g . This critical temperature can be defined 
operationally by the divergence of the relaxation time 
for large systems, or by invoking scaling laws. It is well 
known that at a standard second order transition one 
can define critical exponents which are related to each 
other and which have been shown, by the renormaliza- 
tion group theory, to obey universality rules. These state 
essentially that in a given dimension, all systems sharing 
a few basic parameters (such as the number of order pa- 
rameter components) will have identical values of critical 
exponents. Physically this remarkable regularity arises 
because averaging over larger and larger volumes as the 
critical temperature is approached means that all micro- 



scopic distinctions between systems are blurred out. 

There have been a number of numerical estimates of 
critical exponents in spin glasses. A first point which is 
clear is that the estimated exponents at a spin glass tran- 
sition are very different from those at a standard second 
order transition. For instance, in dimension 3 the spe- 
cific heat exponent a for an Ising ferromagnet is close 
to zero; in the Ising spin glass it is close to —2 J|,|]]. It 
is then important to check if the universality rules hold 
at spin glass transitions as they do at canonical second 
order transitions. This turns out to be difficult to do 
in practice, mainly because the relaxation as the tran- 
sition is approached becomes exceedingly slow, meaning 
that long and tedious numerical anneals must be used to 
make sure that the system is in equilibrium. As a result, 
simulation estimates of exponents have had substantial 
error bars associated with them. 

We have used a combination of scaling rules which al- 
low us in favourable cases to establish accurate values 
of the ordering temperatures together with the static 
exponent rj and the dynamic exponent z (^J^|. Among 
these techniques is that of the scaling form of the non- 
equilibrium relaxation in going from complete disorder 
towards equilibrium, after a numerical quench to the crit- 
ical temperature. This has the advantage of needing no 
preparatory anneal. It would be much more elegant if 
it were possible to establish the ordering temperature 
uniquely from such non- equilibrium measurements, as 
has been done in regular magnetic systems ||. Unfor- 
tunately up to now we have found no simple method of 
doing this. 

The exponent values that we have estimated are more 
precise than but in good agreement with other simulation 
estimates and with exponent values obtained from 

series methods I|. The results clearly indicate that the 
exponents r\ and z for spin glasses in a given dimension 
change with the form of the random interactions. For 
instance, the exponents found with the binomial near 
neighbour interaction distribution ±J are not the same 
as those found with a Gaussian near neighbour interac- 
tion distribution. According to the universality rules, 
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such differences in the interactions should not be perti- 
nent. One would expect, with a change in distribution, 
a change in the critical temperature but not a change in 
the exponents. We can note that in experimental spin 
glass systems the estimates for 77 vary considerably from 
one material to another, again suggesting a breakdown 
of universality |l0| . It seems that ordering at a spin glass 
transition cannot be simply classified with standard sec- 
ond order transitions. This is an important point, as the 
general renormalization argument would seem to be very 
robust so the breakdown is unexpected and goes against 
the conventional wisdom. The empirically observed non- 
universality implies that the physics of critical behaviour 
in spin glasses and other complex systems is much richer 
than in regular systems. It would be very useful to find 
a theoretical underpinning for this result, and to under- 
stand how the exponents change with the various control 
parameters such as the interaction distribution. 

As we have noted above, relaxation at the approach 
to a spin glass transition becomes extremely slow. It 
also becomes highly non-exponential. In a massive sim- 
ulation Ogielski jjj measured accurately the relaxation 
of the autocorrelation function q(t) —< Si(t).Si(0) > 
for large samples of the ±J Ising spin glass in dimen- 
sion 3. At the approach to a continuous transition 
one expects q(t) to take the form t~ x f(t/r) where x 
is a combination of critical exponents and r is a char- 
acteristic time which diverges at the critical tempera- 
ture. Ogielski parametrised his results very succesfully 
using for the scaling function / the stretched exponen- 
tial form, exp[— (t/r) 13 }, with a temperature dependent 
stretched exponential (or Kohlrausch (TTJ ) exponent (3. 
The stretched exponential is used ubiquitously by the 
glass community for fitting experimental relaxation data. 
Having /3 less than 1 can be interpreted as an indication 
of a distribution of relaxation times in the system. In 
Ogielski's data, as the temperature approaches T g and r 
diverges, (3 tends to a value close to 1/3. Numerical data 
on different spin glasses are compatible with similar be- 
haviour in each of the other cases |l2|]l3[| . However good 
the data, it is of course always necessary to extrapolate 
to estimate the limiting value of /3 as r diverges, which 
leads to a margin of uncertainty. Nevertheles the general 
tendency appears to be very similar in all the spin glasses 
which have been studied. 

In experimental spin glass systems it is difficult to 
study the form of the relaxation for a range of temper- 
atures above T g directly because of time scale problems, 
but muons can be used as local probes which are sensi- 
tive to the relaxation of the magnetic spins. The form of 
the decay of the muon polarization with time gives indi- 
rect information on the form of the spin relaxation, and 
it turns out that the muon time window is appropriate. 
Measurements have shown that in a standard metallic 
spin glass the spin relaxation is homogeneous and close 
to exponential only at temperatures of the order of 4 



times T g , and that the muon polarization decay pattern 
is compatible with an Ogielski form of spin relaxation 
for temperatures close to T g jL4],[l5). This means that 
the behaviour observed numerically in Ising spin glass 
simulations is mirrored in experimental, Heisenberg, spin 
glasses. 

This pattern of behaviour is in fact even more general 
and is not confined to spin glasses alone. The relaxation 
in most structural glasses has a characteristic time which 
appears to diverge at a temperature To somewhat below 
the "glass temperature" T g . (In the structural glass con- 
text the latter is conventinally defined as a temperature 
where the relaxation time scale is of the order of a few 
seconds. This is a convenient parameter for defining the 
practical onset of glassy as against liquid behaviour, but 
this glass T g is obviously not a thermodynamic temper- 
ature. To, the divergence temperature, is the equivalent 
of the spin glass freezing temperature). The relaxation 
data are well fitted by stretched exponentials, and for a 
large number of glasses (3{T) drops regularly as T is ap- 
proached, with (3 tending to a value in the region close to 
1/3 when the points are extrapolated to the temperature 
where the relaxation time becomes infinite |l6|-|l8|] . Once 
again, extrapolation is obligatory, as the time needed to 
establish equilibrium before the start of the relaxation 
measurement excedes laboratory time scales well before 
To is attained. Nevertheless there is a striking similarity 
between the behaviour of the relaxation in the Ising spin 
glasses and in these structural glasses, two families of sys- 
tems which have nothing in common on the microscopic 
level. This suggests that in these systems the special 
pattern of the dynamics as the freezing temperature is 
approached is a necessary consequence of the physics of 
the glassy transition. 

We have proposed a possible interpretation couched in 
terms of the morphology of phase space at temperatures 
above the glass transition []l9|j20[ . Let us concentrate on 
the Ising systems. The total phase space of an N-spin 
Ising system consists of 2 N configurations which can be 
mapped onto the 2 N corners of a hypercube of dimension 
N. At finite temperature T some of these configurations 
have such high energies that they are thermodynami- 
cally " inaccessible" . The system performs a random walk 
among the low energy configurations that remain " acces- 
sible" . Relaxation through successive single spin flips is 
simply a random walk where each step takes the system 
from one configuration to a near neighbour configuration 
on the hypercube. The difficult problem is to establish 
the morphology of the set of accessible configurations, 
given the Hamiltonian for the system. Suppose that in 
a glassy system we take an extreme limit and assume 
that a temperature T corresponds to a fraction p of ac- 
cessible sites, distributed at random over the hypercube. 
The relaxation problem is now mapped onto a geomet- 
rical problem. As p drops we will have the closed space 
equivalent of a percolation transition, with a single giant 
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cluster of sites existing down to a critical concentration 
p c Relaxation becomes a random walk among the 

sites of this giant cluster. We have argued f0| , using the 
known behaviour of random walks on a percolation clus- 
ter in Euclidean space |2^ ], that for the high dimension 
hypercube relaxation will approach the stretched expo- 
nential form exp[— [t/r) 1 /$)} as p c is approached. Nu- 
merical simulations on this geometrical model confirm 
this conjecture [l9| , p0| |. 

Our explanation of the particular form of the relax- 
ation as r diverges at a glass transition is then the fol- 
lowing. At the transition the accessible phase space splin- 
ters into many topologically disconnected and inequiva- 
lent pieces, like the small clusters below the critical con- 
centration in a conventional percolation transition. At a 
temperature just above the critical temperature the ac- 
cessible phase space is still topologically connected, but 
it has become so sparse and complicated (fractal-like) 
that the system takes a macroscopic time to explore all 
the nooks and crannies of this structure, in order to be 
able to make up its mind whether the attainable phase 
space is indeed topologically connected or not. As the 
random walk path is tortuous, the relaxation is strongly 
non-exponential, and the hypercube geometry argument 
leads us to a specific form of relaxation. We have dis- 
cussed the Ising phase space explicitly, but any physical 
system has a phase space and the general image should 
carry over to other glassy transitions. This shattering 
of phase space at the transition is very different from 
the description of a second order transition in a regular 
system where at the ordering temperature phase space 
splits up into two (or a small number) of configuration 
clusters related to each other by symmetry. One can sur- 
mise that the non-universality of the critical exponents 
is also a consequence of the geometry of phase space, but 
the connection remains to be established. 

We should at this point make a caveat. The discussion 
as it stands relates to glasses which are known as "frag- 
ile" glasses in the glass terminology, i.e. to glasses where 
To is relatively close to T g . In "strong" glasses, where the 
relaxation time grows following an Arrhenius law as the 
temperature is reduced to T g , it appears to be the ele- 
mentary step time which has increased to a macroscopic 
value at T g , and it is impossible to probe close to Tq 
(supposing it exists) in these glasses. f3 remains close to 
1 down to T g for these cases, which is entirely consistent 
with the general arguments we have given. In intermedi- 
ate cases the extrapolation down to To can only be made 
with a considerable degree of uncertainty. Relating the 
degree of fragility to the molecular structure of each glass 
remains a fundamental problem. 

In conclusion, we have discussed numerical and experi- 
mental data on spin glasses and glasses. The results indi- 
cate empirically that the conventional universality rules 
are not obeyed at spin glass transitions. On the other 
hand, as the critical temperature (where the relaxation 



time scale diverges) is approached, the form of the re- 
laxation in spin glasses or fragile structural glasses tends 
towards a stretched exponential with a Kohlrausch ex- 
ponent close to 1/3. We have provided a general phase 
space image that could explain these observations. 
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FIG. 1. The dynamic exponent z and the critical exponent 
r\ as functions of the critical temperature T g for four Ising spin 
glasses in dimension 4. The sets of interactions, from left to 
right, are exponential, gaussian, uniform and binomial M. 
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FIG. 3. The Kohlrausch exponent /3(T) as a function 
of (log(r(T) lO 1 ^})" 1 obtained from dielectric relaxation 
measurements mjt on a selection of polymers (PSF (■), 
PViBE(v), PE (0), PEI (T), PH (♦) and PMA (□)) and 
on the 3D Ising spin glass model with binomial interaction 
distribution (+) Q. The dashed line represents the practical 
limit (r ~ 4 months) for equilibrium measurements. All the 
sets of data are compatible with an extrapolated value of j3 
near 1/3 as r goes to infinity (the lines through j3 — 1/3 are 
to guide the eye). 




1.4 1.6 1.8 2.0 

Temperature 



FIG. 2. The Kohlrausch exponent j3 for the relaxation of 
the autocorrelation function as a function of T, for the 3D 
Ising spin glass model with binomial interaction distribution 
M. The points extrapolate to close to 1/3 at the critical 
temperature T g ~ 1.175J (dashed line). 
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